Directed flow in non-adiabatic stochastic pumps 
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We analyze the operation of a molecular machine driven by the non-adiabatic variation of external 
parameters. We derive a formula for the integrated flow from one configuration to another, obtain 
a "no-pumping theorem" for cyclic processes with thermally activated transitions, and show that in 
the adiabatic limit the pumped current is given by a geometric expression. 
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Assemblies of molecules that perform specific tasks, 
known as molecular machines or motors, are ubiquitous 
in biological systems they act as pumps across cell 
membranes, carry loads in cells and cause muscles to con- 
tract. Remarkably, the construction and manipulation 
of artificial molecular machines has in recent years be- 
come feasible [2J. Theoretical models, such as thermal 
ratchets, have provided an understanding of the micro- 
scopic mechanisms underlying the operation of molecular 
machines Q. These systems and models all exhibit the 
rather surprising ability to produce directed motion or 
current in the face of the violent thermal agitation of the 
mesoscopic world [4J. 

Stochastic pumps are molecular machines in which di- 
rected current is produced by varying external param- 
eters, such as chemical environment or applied fields. 
These come in two varieties: open pumps are active con- 
duits for particle flow between reservoirs, while closed 
pumps generate internal currents, such as directional ro- 
tation in catenanes, mechanically interlinked ring- like 
molecules @. Successful theories have been developed 
to describe stochastic pumps under adiabatic condi- 
tions, that is when the external parameters are driven 
slowly B, 0, @, HI, or for weak oscillatory perturba- 
tions [10( 1 . However, a general theory of non-adiabatic 
pumps is lacking. In this Letter we formulate and ana- 
lyze a generic model of a non-adiabatic, closed pump, de- 
scribed by transitions among a network of N states (e.g. 
molecular conformations). We derive an expression for 
the integrated, directed current along an arbitrary link 
in this network (Eq. [4|) and explore its consequences. In 
particular, we obtain a "no-pumping theorem" that ap- 
plies to the cyclic variation of thermally activated tran- 
sitions. 

We will consider a system with a set of configurations, 
or states, labeled i = 1, • • • , N, and will model transitions 
among these states as a Markov jump process. Letting 
Pi(t) denote the probability to be in state i at time t, the 
system obeys the master equation (llj 

p = TZp, (1) 

where p = (pi,-- - ,Pn) and TZ is a transition matrix: 



when the system is in state j, the probability rate to 



jump to state i is R^ > 0; and Rj 



Liil; 

At time t, the probability current from state j to state i 



Jij(t) = RijPj - RjiPi 



(2) 



We imagine that the transition rates depend on a vector 
of external parameters A under our direct control, i.e. 
TZ = 7Z(\). For simplicity, we further assume that Rij = 
if and only if Rji — 0; and that for any A, there is a 
unique stationary state p s (A) satisfying TZp s ~ 0, with 

stationary currents J^-(A) = RijPj — RjiPi- If Jfjfi) = 
for all i,j, then the dynamics are said to satisfy detailed 
balance, and we interpret the stationary distribution to 
be the equilibrium distribution: p s (A) = p e? (A). 

It will prove convenient to define a branching matrix Q, 
obtained by rescaling each column of TZ: Qij — Rij / \ Rjj \ . 
The diagonal elements of Q are all —1, and the off- 
diagonal elements are branching fractions: when the sys- 
tem makes a transition out of state j, then Qij is the 
probability that the transition is into state i. 

As an illustrative example, for which the dynamics sat- 
isfy detailed balance, imagine thermally activated tran- 
sitions between potential wells separated by energetic 
barriers, with well depths Ej(X) and barrier energies 
Bij{\) (Fig. [[J. We then have p* = pf cx e' ^ < and 
transition rates take the familiar Arrhenius form [lij 
Rij — kexp[—/3(Bij — Ej)], from which it follows that 
the branching fractions depend on the barrier energies 
but not on the well depths: 



Qi 



(3) 



We will use this result in our later analysis. 

We will consider a process during which the system 
evolves as the parameters are varied externally, A = A t , 
from Ao = A to A r = B. For a fixed but arbitrary pair of 
states to, n, the integrated current, $ m „ = J Q dtJ mn (t), 
represents the net transfer of probability from n to to. 
This is a measure of the directed motion produced along 



2 



this particular link in the network of states. We will 
show that the integrated current is given by the compact 
expression, 



4>, 



dt 



J^ n (X t ) + V mn (X t )-p(t) 



(4) 



where the vector field V m „(A) = (y mn ,i, ■ ■ ■ ,V mn ,N) is 
given by Eq. [TO] below. Since in general we must inte- 
grate Eq. [1] to obtain p(i), Eq. [4] is not a shortcut for 
calculating $ mn . Rather, it provides a useful and non- 
trivial decomposition of $ into a stationary contribution, 
$ s = J J s dt, due to a baseline current that flows even 
at fixed A; and the remaining excess, or pumped con- 
tribution, & ex , due to the redistribution of probabilities 
induced by the variation of A. The latter is the contrac- 
tion of a path p(t) along the field V mn (A): 



V mn (A t ) • dp(t). 



(5) 



While Eq. [4] is valid quite generally, we note two conse- 
quences that follow in particular situations. 

(1) When A is varied slowly, the system remains near 
the stationary state, p(t) w p s (A t ). This suggests we 
replace dp in Eq.[5]by dAoVp s = £ M dX^dp 3 /dX^) [lj], 
in the adiabatic limit. With this replacement, which can 
be justified by appeal to an adiabatic perturbation theory 
[3], we get 



A mn (X) o d\, 



(0) 



where A mn (X) = V mn • Vp s . This expression is geomet- 
ric: time no longer appears here, and is simply a line 
integral of A mn along a path in parameter-space. Similar 
geometric results have been obtained by Astumian for a 
three-state system Q , and by Sinitsyn and Nememan Q 
for open stochastic pumps. 

(2) When ft describes thermally activated (Arrhenius) 
transitions over barriers, with parameter-dependent well 
depths and barrier energies (Fig. [l}, then we will show 
that Eq. [4] leads to a surprising "no-pumping theorem" 
for cyclic processes: § ex = if the well depths are varied 
while the barrier energies are held fixed, or (trivially) 
vice-versa. Only by varying both well depths and barrier 
energies during a pumping cycle, can we generate a net 
transfer of probability between states. 

Eq. |4] is derived by eliminating p from Eqs. Q] and [2] 
to obtain J mn {t) — J^ m + Y mn ■ p. We now sketch 
the steps of this derivation, skipping tedious but rou- 
tine linear- algebraic manipulations. The final results give 
V m „ in terms of readily evaluated minors of the matri- 
ces ft (Eq. [10]) or Q (Eq. I lib]) . We note before pro- 
ceeding that there is a gauge freedom at play here: since 
J2j Pj — by conservation of probability, Eq. 2] is unaf- 



1 = (1, 1, ■ • ■ , 1) and / is an arbitrary function. Thus 
our results for V mrl are not unique, merely convenient. 

Eq. [lj is a set of linear equations, which we label 
ei, ••■ , e^v • Since detlZ = 0, these are linearly depen- 
dent (one of them is redundant) and ft cannot be in- 
verted to solve for p in terms of p. Specifically, for a 
given p, if p satisfies Eq. [JJthen so does p( Q ) = p + ap s , 
for any value of a. We remove this degeneracy by impos- 
ing the normalization condition 1 • p = 1 , which we label 
e': replacing ejv by e' in Eq. [TJ we get a set of linearly 
independent equations 



P' = ^'p, 



(7) 



where p' = (pi, • • • ,pn-i, 1), and ft' is obtained by sub- 
stituting the vector 1 for the iV'th row of ft. Since 



det TZ' 7^ 0, we solve for p using Cramer's rule [151 ]: 



Pj = det ft'/ det ft', 



(8) 



where TZ'j is obtained from TZ' by replacing the j'th col- 
umn by p'. Expanding det TZ'j along this column (for 
j = m, n), then substituting Eq. [8] into Eq. O we get 



N 

E 

k=l 



(-i) fe K 



(9) 



Here cy, fc (A) = (~l) j Rij r'^k; j)/r', where r' = det ft', 
and r^(o;6) denotes the (a; b) minor of ftj, that is the 
determinant of the matrix obtained by deleting row a and 
column b of TZ' i . Comparing Eq. [9] with the integrand in 
Eq. [H and recognizing that J mn = J^ n when p = p s 
(i.e. when p = 0), we obtain 



Vmn.k — (— 1)*"' (1 — &kN) (o~mn,k 



(10) 



which gives V mn (A) in terms of the elements of ft (A). 

Let us now separately analyze the case in which the 
dynamics satisfy detailed balance, for all A. Recalling 
that Eq. [TJ supports a family of solutions, T = {p^ — 
p + ap s }, we note that if we formally replace p by p( Q ) in 
Eq.[2l we obtain jf"' = J i3 -+aJ^. In general, is not 
a physically meaningful quantity. However, if ft satisfies 
detailed balance, as we assume in this paragraph, then 
Jf = 0, and therefore all solutions p( Q ) G T give the 



same current, J, 



(a) 



Jij , upon substitution into Eq. [2] 



fected by the replacement V ?i 



/ (A)l, where 



We now exploit this observation: rather than solving for 
the vector pgf that satisfies normalization (e'), as done 
above, we instead choose the vector pg J that satisfies 
Pm = (we label this condition e") and then we sub- 
stitute p into Eq. [2] to determine J mn . This approach 
leads to an expression for V m „ that (unlike Eq. [TO]) de- 
pends only on the branching fractions Q and not on the 
transition rates ft. In detail, to solve for p we break 
the degeneracy of Eq. [JJby replacing i m with e", rather 
than e^v with e' as earlier. In lieu of Eq. [7] we now have 
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p" = 11" p, where p" = (pi, • • ■ ,p m - 1 ,0,p m+1 , ■ ■ ■ ,p N ), 
and 1Z" is denned by replacing the m'th row of 1Z with 
(0 ■ • • 1 • ■ ■ 0), where the 1 is on the diagonal. Now, 
(1) using Cramer's rule to solve for p, then (2) taking 
J mn = RmnPti (since p m = 0), and (3) recalling that 
Qij = Rij/\Rjj\, after some effort we obtain 



dtV mn (X t ) ■ p(i) 



(11a) 



w = ( -irA ( i-u"^iib) 

q{m;m) 

Here q(m; m) is the (m; m) minor of Q, and similarly 
g(m, k; m, n) is the determinant obtained after deleting 
rows m and k and columns m and n of Q; finally A = — 1 
if n < m < k or k < m < n, otherwise A = +1. 

Now let us use Eq. [11] to establish a "no-pumping the- 
orem" pertaining to cyclic processes. We continue to as- 
sume that TZ(X) satisfies detailed balance for all A, and we 
picture the dynamics as arising from thermally activated 
transitions over barriers, with externally controlled well 
depths and barrier energies (Fig. [l}. We now imagine a 
process during which the parameters are held fixed at A ,4 
prior to t = 0; then during the interval < t < T they are 
made to trace out a closed loop in parameter space, after 
which they are again held fixed at Xa- In this scenario, 
from the distant past to the distant future the vector 
p(t) also evolves through a closed path, returning to its 
initial equilibrium state: p(±oo) = p eq (\A)- We are in- 
terested in the integrated current, & m n = J^^dtJmn 
during such a process. Let us separately consider two 
cases. First, if only the barrier energies (and not the 
well depths) are varied, then the system simply remains 
in equilibrium at all times, p(t) = p C9 (Aa), and there 
are no currents. Now consider the less obvious case in 
which the barrier energies are held fixed, but non-zero 
currents J mn (t) are generated through a cyclic varia- 
tion of the well depths. Since V„ m is determined by 
the elements of Q (Eq. Illb[) . which in turn depend only 
on the fixed barrier energies (Eq. [3]) , Eq. Illal becomes 
^ mn = V m „ • dtp(t) = 0. Thus, in order to gener- 
ate non-zero integrated current during a cyclic process, 
both well depths and barrier energies must be varied. 

This result extends to encompass repeated, periodic 
cycling of the parameters, A(t + r) = A(t), in which case 
Floquet theory ensures that the system relaxes to a time- 
periodic state, p(t + t) = p(t) 16]. The arguments of 
the previous paragraph then apply to this periodic state: 
if the barrier energies are held fixed, then V mn (At) is 
constant in time, and therefore Eq. Illal gives $ mn = 
over one period of driving. 

We have just argued that all $ mn 's vanish in a cyclic 
process during which Q is fixed and detailed balance is 
satisfied. It is of interest to count the number of con- 
straints represented by these conditions, and to compare 
this with the general number of constraints needed to en- 



sure absence of directed flow. Assuming E links among 
our network of N states, there are 2E non-zero branch- 
ing fractions Qij. Such a network can be decomposed 
into E — N + 1 closed loops, with each link belonging 
to at least one loop 17 1. The 2E non-zero branching 
fractions are not independent: conservation of probabil- 
ity imposes N conditions YliQij = (j = 1, • • • ,N), 
and detailed balance imposes an additional E — N + 1 
conditions, representing the absence of thermodynamic 
force around each closed loop [3] ■ Therefore by choosing 
E — 1 specific branching fractions to be A-independcnt, 
we guarantee that the others will also be A-independent. 
On the other hand, we can ensure that all <S> m nS vanish 
by insisting that there be zero integrated current around 
every closed loop in the network decomposition. Since 
there are E — N + 1 such loops, it is clear that fixing 
the branching fractions (i.e. imposing E — 1 constraints) 
is not the most general condition for the absence of in- 
tegrated flow. However, it is a simple condition, which 
has the same leading order scaling for systems with many 
states, assuming E oc N 2 . 

We now illustrate our results using a model system, 
motivated by an experiment by Leigh et. al. J5( and 
analyzed in the adiabatic limit by Astumian [9]. We 
consider thermally activated transitions among three 
states, depicted by the wells in Fig. [lj with rates Rij = 
feexp[— /3(Bij — Ej)], and we will take k, /3 = 1 to set 
the units of time and energy. TZ satisfies detailed bal- 
ance, but by varying the well depths and barrier energies 
we can induce non-zero currents. Recalling Eq. [3l and 
defining ipx = exp(-5i 2 - B 13 ), i/j 2 = exp(-B 12 - B 23 ), 
tp3 = cxp(— £?i3 — -B23), and K = J^j^Pji we evaluate 
Eq. Illbl for (m,n) = (2, 1) to obtain 

V 21 =K- 1 (-^ 1 -^ 2 , 0,-Vi) ^ K-\-^ 2 ,^ u 0), (12) 

where in the last step we have used the gauge freedom 
Vai V21 + (fa/K) 1. 

Energy 




Configuration 

FIG. 1: A model stochastic pump satisfying detailed balance. 
The particle makes thermal transitions among three states 
with energies Ej(X), over barriers with energies Bij(X). These 
are varied with time to induce currents. 

When A is varied adiabatically around a closed path, 
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the pumped current is given by Eq. [6j with A21 = 
{—ip2^p e i + i'i^P2 i )lK. If the barrier energies are held 
fixed during this process, then ipi, "0 2 > and K are con- 
stant, hence the integrand is a total differential and there 
is no pumped current, as predicted in Ref. |9|. 

Now let us analyze cyclic but non-adiabatic variation 
of the well depths and barrier energies. We first con- 
sider a process during which the barriers are held fixed. 
Specifically, we take (-B 12 , B 2 3, -B13) = (—0.3,0.5,0), and 



E {t) = -2- 



cos 



27T 1 - 



j-1 



(13) 



for < t < T = 10. Thus the well depths E 3 (t) undergo 
one cycle of pumping, with phases staggered by 27r/3 in 
a piston-like sequence. Outside this time interval all pa- 
rameters are fixed, so the system ultimately relaxes to its 
initial equilibrium state. The solid line in Fig. [2] shows 
the integrated current $21 (t) = Jq dt J%x(t) , obtained 
by numerical integration of Eqs. Q] and O We see that 
probability sloshes back and forth on the link between 
states 1 and 2: initially there is a gentle flow from 1 to 
2 (g?<I>2i/g?t > for r < 2), then an interval of stronger 
current in the opposite direction, followed by another re- 
versal shortly before r = 7.5. The eventual decay of <5> 2 i 
to zero indicates a net cancellation of these flows, as pre- 
dicted by our no-pumping result. We next consider a pro- 
cess during which both well depths and barrier energies 
are varied with time: the E^s are again driven according 
to Eq. 1131 but now each barrier moves in synchrony with 
the well to its immediate right in Fig. [TJ e.g. as E\ goes 
down and then up, so does -B13, so that their difference 
remains fixed at B13 — E\ = 2 = B\ 2 — E% = B23 — £3. 
The integrated current $21 (t) is shown by the dashed 
line in Fig. [2j the asymptotic value $21 ~ 0.1 reveals 
a net transfer of probability from state 1 to state 2 over 
the cycle. Note that in both cases non- vanishing currents 
persist for some time after r = T, reflecting the decay to 
equilibrium that occurs after the parameters stop being 
varied. 

Our no-pumping theorem applies to a single particle 
jumping among potential wells (Fig. [I]). When more par- 
ticles are present, and they interact with one another, 
the general results Eqs.mDIB and I llbl remain valid, but 
now the roles of wells and barriers are played by many- 
body energies. We then find non-zero cyclic currents 
even when the single particle barriers are held fixed [3] • 
These results are consistent with the experimental ob- 
servation of currents in 3- but not 2-catenanes H; these 
linked-ring molecules are naturally modeled as two- or 
one-particle systems, respectively, with the particle(s) 
jumping among binding sites whose affinities are varied 
externally [9j. This agreement suggests that the non- 
adiabatic framework described in this Letter will prove 
useful in the design and analysis of nanoscale stochastic 
pumps. 
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FIG. 2: The integrated current $21 for non-adiabatic cycles 
with fixed barriers (solid line) or varying barriers (dashed). 



Recently, we have learned of a generalized no-pumping 
theorem derived by Chernyak and Sinitsyn [191 ]. Their 
results apply to systems satisfying detailed balance and 
account for the topology of the network. 

We gratefully acknowledge Doran Bennett and Nikolai 
Sinitsyn for useful discussion and correspondence, and 
the University of Maryland for financial support. 
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